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We construct a topological invariant of the renormalization group trajec- 
tories of a large class of 2D quantum integrable models, described by the 
iSJ ' thermodynamic Bethe ansatz approach. A geometrical description of this 

fvj , invariant in terms of triangulations of three-dimensional manifolds is pro- 

C^ • posed and associated dilogarithm identities are proven. 

a^ 

1 Introduction 

Oh: 

O . Many quantum integrable 2D models can be considered as perturbations of a Con- 

formal Field Theory. The RG evolution of these models is described, in the thermo- 
dynamic Bethe ansatz (TBA) approach [^ by the ground state energy E{R) of the 
^ ■ system in a infinite cylinder of radius R, through an equation of the form E{R) = 

■^ ^^^^ /(i6't'a(^) log (1 -|- ya(6')), where 9 is the rapidity, a labels the pseudoener- 
gies (cf. Q), Ya{6) are i?-dependent functions, determined by a set of coupled integral 
equations known as TBA equations; the f^'s are elementary functions of R and ^, re- 
lated to the asymptotic behaviour of the y^(^)'s for 9 -^ ±00. An important observation 
of Al.B.Zamolodchikov|^ was that any solution {Ya{9)} of the TBA equations satisfies 
a set of simple functional equations, called the Y-system. Conversely, it can be shown 
that under appropriate analytic and asymptotic restrictions the Y-system is equivalent 
to the set TBA equations. For instance, the Y-system describing the flow from the 
tricritical Ising model to the critical Ising one is given by 

77]- 777 fjr 777 

Y^{9 + -)Y^{9--) = 1 + Y2{9), Y2{9 + -)Y2{9--) = 1 + Y,{9) . (1.1) 



This system, like any other Y-system, can be considered as a recursion relation in the 
imaginary rapidity. To simplify the notations we can put T„(i?, 9) = Y- 



then the eq. (LI) becomes 
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1 + T„, . 



(1.2) 



It is immediate to verify that this Y-system is periodic: Tn+5 = T„, i.e. Yi{6 + |i7r) = 
¥2(9) , 12(^+1^71") = ^i(^) • In ref. Q it was also conjectured that a similar periodicity 
Ya{9 + Piir) = Ya{9) where the period P is related to the dimension of the perturbing 
operator by A = 1 — 1/P. In the UV limit i? — > the vacuum energy E[R) is expected 
to behave as E{R) ~ —7rc/6R, where c is the effective central charge of th UV fixed 
point. In general one finds 
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1 + y„(oo) 



(1.3) 



where L{x) is the Rogers dilogarithm. It can be defined as the unique function three 
times differentiable which satisfies the following five term relationship, known also as 
the Abel functional equation (see for instance 0) 



L{x) + L{y) + L ( -!—^] +L{l-xy) + L^ ^ ^ 



1 — xy 



1 — xy 



3L(1) , 



(1.4) 



with < x,y < 1 and the normalization -^(1) = ^. This equation has a ^5 symmetry 
[15, 1^ which becomes manifest rewriting it in the form ||^] 



n=l ^ ' 



3L(1) , 



(1.5) 



where the T„'s fulfil the recursion relation (1.2). Clearly eq. (1.5) defines an invariant of 
the RG flow of the tricritical Ising system, being the L arguments non-trivial functions 
of R and 0, while the r.h.s. is constant along the RG trajectory. This behaviour suggests 
[0] the existence of large families of functional equations, which are generalizations of 
( |1.5D . The simplest examples are classified in terms of an ordered pair G x H oi ADET 
(T„ = A2n/'^2) Dynkin diagrams M 
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Ti{n - l)T^(n + 1) = J{{1 + T^(n))^- \{[l + 
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Here Gac and -ff^^ are the adjacency matrices of the appropriate Dynkin diagrams 
whereas re and rn are the ranks of the corresponding algebras. The solutions to this 
system are conjectured to be periodic functions: T^(n + P) = T|.(n) , P = h + g , with 



a and b the conjugate nodes to a and b, while h and g are the dual Coxeter numbers. 
Then, the T's solving eq. ( |1.6| ) are, according to our conjecture, arguments of functional 
equations for the Rogers dilogarithm, L{x), of the following form 

a=l 6=1 n=0 ^ a\ ) / 

Notice that since L(x) is a multi-valued function, the r.h.s of eq. ( |1.7D cannot be fixed 
unambiguously, unless a consistent choice of the sheet of the associated Riemann surface 
is implemented. The quoted value corresponds to solutions of the recursion relations 
with all the T's > 0. It is interesting that the Ajq x A\ {N = 1,2,3) identities re- 
duce respectively to the Euler, Abel and Newman functional equations (cf. |^] ) when 
symmetric conditions on the conjugated nodes are chosen. 

Other examples of Y-systems with their corresponding dilogarithm identities were 
found and discussed in ref. Q. Furthermore, two new families of Y-systems and dilog- 
arithm functional equations, classified in terms of a pair of co-prime integer (p, q) asso- 
ciated to the sine-Gordon model at rational points and to its reductions, were proposed 
and partially studied [P, |lO| (see [|l^] for even more exotic relations). Very soon af- 
ter their discovery, the general solution of the Aj<; x A\ systems was found and the 
corresponding dilogarithm identities independently proven in [10, 13 1. We report here 



the solution proposed in ref. ||lO|, which was related to the fact, already known to 
Lobachevskij , that dilogarithms define volumes of particular tetrahedra in hyperbolic 
space and other manifolds. The seed of the idea was introduced in 1992 in the field of 
integrable models by ref. |^. Other ideas were then expound and elaborated in ref. [p!^], 
where a manifold on which Rogers dilogarithms define volumes of ideal tetrahedra was 
identified. 

2 The geometrical picture 

The starting point of ref. [1^] was the assumption that the Abel equation mimics in 
some way a volume calculation of a compact manifold. Although the existence and the 
nature of such a manifold could not be further argued, this geometrical point of view 
was rather powerful and allowed to prove in a rigorous way, just using eq. ( |1.4[) as the 
only ingredient, the Aj<s x A\ identities and to solve the Y system of a general class 
of thermally perturbed minimal models. Let us mention that after the main results of 
ref.s [0, [l^, 1^] were already appeared, we became aware Q of an undated manuscript by 
D.Zagier where this geometric approach is clarified and simplified. We can now rephrase 



* We are grateful to the referee of an earlier version of fid] for a copy of Zagier's notes. 



our construction as follows. Consider an arbitrary triangulation of an arbitrary three- 
dimensional manifold, made by M oriented tetrahedra {Tj}. The fact that this set forms 
a triangulation means that the boundary of the associated 3-chain is zero: 



M 



diJ2Ti) = 0, (2i 



where the boundary of the tetrahedron T = [ABCD] is defined in the standard way: 

d [ABCD] = [ABC] - [BCD] + [CDA] - [DAB] . (2.9) 

Now associate to each vertex i [i = 1, . . . N) a real number f\ Xj. Through this map we 
can assign to each tetrahedron T, the cross-ratio t = (abed) of the values associated to 
its vertices 

T = [ABCD] -^ t = (abed) , (abed) = i" ~ ''?,1^ ~ "^j . (2.10) 

[a — d){b — c) 

In the following we shall use some elementary properties of the cross-ratios: 

(abed) = (cdab) = (bade) = 1 — (acbd) = (bacd)^^ 

(aabc) = 1 , (abbe) = oo , (abcb) = . (2.11) 



We shall prove that 



M 



Y,L{t,) = nL{l) , (2.12) 



where the integer n is unambiguously fixed when the sequence of real numbers {xj} is 
ordered in such a way that all the ij's belong to the interval [0, 1]. Since the number 
M of tetrahedra is less or equal to the number N of vertices, there are simple algebraic 
relations among the tj's. For suitable triangulations they coincide with the Y-system 
of known integrable models. Before giving a general proof of eq. ( |2.12| ) let us see a 
couple of simple instances. Consider a polytope made of tetrahedra in four dimensional 
Euclidean space. Well known examples of this kind are the three regular polytopes 5-cell 
(or 4-simplex) , 16-cell (dual to hypercube) and 600-cell. They are all triangulations of 
5*3 (see for example pT|]). The 5-cell of fig. la is formed by the following five tetrahedra 



Ti = [BACD], T2 = [CBDE], Tg = [DCEA], T^ = [EDAB], Tg = [AEBC] . 
The pentagonal relation (|1.4|) is then secured assigning five real numbers {a <b < c < 



^More generally, the XiS ma y be chosen as points of an arbitrary circle of the complex plane in order 



to have real cross-ratios in eq. (2.10) 





lb) dual polytope. 



d < e} to the vertex set {A, B, C, D, E}; one immediately end-up with [| 
L{{hacd)) + L{{chde)) + L{{dcea)) + L{{edab)) + L{{aebc)) = 3L(1) . 



(2.13) 



It is instructive to draw the adjacency graph of this set of tetrahedra; associating each 
tetrahedron to a vertex an connecting with a Une two tetrahedra sharing the same face, 
we end up with the graph of the dual polytope, which in this case is again a 5-cell 
(see fig. lb). For each loop of the adjacency graph one can write an algebraic relation 
among the cross-ratios of the tetrahedra touched by the loop. In particular, it is easy to 



verify that the L arguments fulfil the recursion relation (1.2). Moreover the well known 
identity L{x) + L[\ — x) = -L(l), known as Euler equation, can be written as 



L{{abcd)) + L{{acbd)) = L(l) , (2.14) 



hence 



L{{bcad)) + L{{cdbe)) + L{{deca)) + L{{eadb)) + L{{abec)) = 2 L(l) . (2.15) 

In the following we shall use the short-hand notation L{abcde} = 3-^(1) for the identity 
(|2.13|) and L{abcde} = 2L(1) for the identity (|2.15| ). As a more engaging example, let 
us consider the 16-cell. This polytope has 8 vertices, characterized by the property that 
each vertex is linked to all the others but one (see fig. 2a). Ordering the vertices in such 

•'■Because of the projective invariance of the cross-ratio, three of these values can be fixed arbi trarily, 
then these five cross-ratios depend actually only on two free parameters, as required by eq. (|l.4|). 



a way that the unhnked ones are given by the pair i, i + A mod 8, the 16 tetrahedra 
spht in two sets, associated to the foUowing cross-ratios: 



ti 



{xjXj^lXj+lXj+2) , Uj = {xjXj^iXj^2Xj^3), (j = l,2. 



(2.16) 



where all the indices are taken modulo 8 and xi < X2 < . . . xg. Starting from the sum 
of the following five-term relations written as 

22 L(l) = L{xiX2X3X4X5} -I- L{xiX3X4X5Xq} + L{xiX2X4X5X7} + L{xiX2X3X5X8} 
+ L{x5XqX7XsXi} + L{x5X7XsXiX2} + L{x5XqXsXiX3} + Llx^X^X^XiX^} 

one can easily verify that each five-term relation contains two tetrahedra of the trian- 
gulation ( 2.16| ); the remaining 24 tetrahedra combine in pairs according to eq. ( |2.14| ), 
so one is directly led to the sought-after 16-term relation 



Y,[Hk)+Liui)] = 10L{l) 



(2.17) 



i=l 



Drawing the adjacency graph of these tetrahedra we find the dual polytope, i.e. the 
8-cell (or hypercube) represented in fig. 2b. There is an algebraic relation for any set 




2a) 16 cell. 




of tetrahedra connected by a closed path of such a graph. In order to find a link with 
known Y-systems it is convenient to transform the argument t = (abed) of L, using the 
the relation t = —{adch) = t/(l — t). Define now Yi{n) = tsn+s and Y2{n) = u^n, where 



all the indices are understood modulo 8. One can check that the following relations are 
fulfilled 

Yi{n - l)Yi{n + 2) = {1 + Y2in)){l + Y2in + 1)) 

Y2{n-l)Y2in + l) = (l + Yi{n))/il + 1/Y2{n)) . (2.18) 

This is the Y-system corresponding to the (/)i3-thernial perturbation of the Mr,^ min- 
imal model. While the discovery of the Y-system associated to the 600-cell is still an 
open challenge, two infinite families of 5*3 triangulations associated to known integrable 
models may be found in ref. [|l^] . Dilogarithm functions obey a general property known 
as "beta-map" [^, ^, ^, that we will shortly introduce. The beta-map and various 
geometrical notions relating dilogarithms to three-dimensional compact manifolds will 
merge together giving a more unified point of view. Swapping between these two main 
tools we will then be able to prove many of the ADE-related dilogarithm identities. The 
beta-map condition adapted to the Rogers dilogarithm states that a sufficient condition 
for having the dilogarithm functional equations 

J^L(/C,)gZL(1) , (2.19) 

i 

is that the following condition is satisfied 

^/C,A(1-/Ci) = , (2.20) 

i 

where the wedge product is defined through the properties 

xAz = —zAx ,{xy) A z = X A z + y A z . (2-21) 



For a more detailed discussion of the relation ( |2.ig| - 2.20| ) the reader is referred to [p!3| , 0, 



|T5|. If the L argument is expressed as a cross-ratio t, we have t A (1 — t) = (abed) A (acbd) 
and 

{abed) A (acbd) = (abc) - (bed) + {cda) - (dab) , (2.22) 

where we have defined 

(abc) = (a-c) A{a-b) + {a-c) A{c-b) + (b-c) A{a-b) . 



Comparison between eq. (|2.9| ) and eq. (|2.22 ) shows that the beta map acts on the cross- 



ratios associated to the tetrahedra like the boundary operator d, hence we have the 
implication 

dC£Tj) = ^ ^t,A(l-t,) = , (2.23) 



which proves eq. ( 2.12D . It is important to note that the impUcation ( 2.23| ) cannot 



be inverted. Actually we found Y-systems which do not correspond strictly to a tri- 
angulation. The simplest example is the solution of the D^ = D^ x Ai systems of 
eq. (P): 



N 

T,(n + l)T„(n-l) = J];(l + T,(n))^i^^ . (2.24) 

6=1 

These systems emerge from the TBA analysis of the sine-Gordon model at the reflection- 
less points ^ = 1/(A^ — 1) in the attractive regime, or equivalently (with the change 
T — > 1/T) at the points ^ = {N — 1) [IC] in the repulsive regime. The solution will 



be written in terms of cross-ratios of a unique function j ^ hj (j = 1, . . . N) exactly 
as in the models treated previously, but now h is quasi-periodic, in the sense that 
hj-N = qhj +p, where p and q are real numbers. For the ^Ar-2 tail of the diagram (see 
fig. 3a) we have 

Tain) = -{hjh,hi+ihj+i) , i = [^ttl] ^ J- = [n±a] , a<N-2 , (2.25) 

with the square brackets indicating that the indices must be taken modulo N. For the 
two nodes A^ — 1 and A^ (the "fork" of the diagram) we can write 

Tiv(n±l) = Z±(n±l) , Tiv-i(n±l) =Z^(n±l) , (2.26) 

with Z^{n) = — (/ij/ij/ij_|_ioo), Z^{n) = — (/ij+i/ij-|_i/ijOo), where i = [ "~2~ ] and j = 
[ ^~^ 2 ~ ]• ^^ ^^ ^°^ easy to verify the conjectured periodicity Ta(n + 2A^) = Ta{n). 
Introducing, besides the A^ points associated to hj, N further vertices j = 1,...N 
associated to the values j -^ qhj +p and a vertex associated to cxd, it is possible to 
represent this Y-system with a 3-chain D]^ of A^^ tetrahedra and 2A^ + 1 vertices. Now 
the boundary of this chain does not vanish, however we have d{D]\f + B^} = , where 
Bjsf is the auxiliary 3-chain 

N-2 

We can then write a dilog identity associated to the triangulation D^ + B^ (which 
turns out to be again a triangulation of 53). Since we have 

(hjhNhj+ioo) + (hjhj^hjjoo) = 1 , 

the tetrahedra of Bj\f can be paired off using the Euler identity L(x) + L(l — x) = L{1), 
so we end up with the result conjectured in eq. (|1.7| ). Note that because of the pairing 
just mentioned, the chain B]\f gives a vanishing contribution to the beta map ( |2.20| ), 
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N-2 
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3a) Dm diagram. 



Y,(4) 




Y,(0) 



3b) Adjacency graph of D4. 



hence also the chain D^Vi although it is not a triangulation, fulfills this condition. As 
an explicit example, the adjacency graph of the D4 system is drawn in fig. 3b; the grey 
circles denote the tetrahedra of the B4 chain and the dotted lines correspond to their 
faces. Other infinite families of such "quasi-triangulations" can be found in ref. ||l^ 



3 General ADE x ADE identities 

In this section we will generalize the proof of the functional dilogarithm identities to 
more general cases, where the existence of an associated triangulation is not known, 
using the only assumption that the associated Y system is periodic. Let's introduce an 
auxiliary set of equations again classified in terms of a pair of Dynkin diagrams (G, H) 



TH 



re 



T^in + l)Tl{n - 1) = J] T^i^)""'" + 11^^ 



•bi^-xGac 



(3.28) 



d=l 



c=l 



When specialized to Aj^ x ADE the three term relations ( 3.28| ) coincide with the well 
known (cf. |1^]) fusion-relations appearing in integrable lattice models. They are often 
called T— systems. Analytically for the lower rank cases and numerically more generally 
one can again show that they fulfil the cyclic property 



T',{n + P)=T^{n) , P = h + g . 



(3.29) 



Furthermore defining 

re TH 



Ya'in) = n T'Anf- JJ T^{ny^'^ , (3.30) 

c=l d=l 



and 



1 + Y^{n) = Tl{n + l)Tl{n - 1) J] ^.^ W"^- , (3.31) 



d=l 



one can easily show that the functions Y defined in (3.30) satisfy eq.s ( |1.6| ). Now 



writing the y's in terms of the T's using ( 3.30| ) and the periodicity ( 3.29| ) one can show 
that ( p. 20 ) is vanishing. We have 



P rjj ro 

n=l a=l 6=1 

ra rn P r„ 

{n,a,fe} c=l <i=l m=l e=l 

where for convenience the incidence matrix A^ ' of an Affine Ap-i Dynkin diagram has 
been introduced. Using abf\cd = af\c + h/\c + af\d + h/\d and 

dim{M) dim{K) dim{M) dim{K) 

E E ( n ^cY"^)^( n ^-") = E A4,i^e5(Xe6AX,,) = 
a=l 6 c=l e=l {a,6,c,e} 

= I J^ MacKe6(^c6 A Xae " X^e A X^fe) = , (3.33) 

{a,b,c,e} 



for M and i^ arbitrary symmetric matrices we conclude that eq. (3.32) vanishes and 
therefore we have dilogarithm identities. 

4 Conclusions 

We have described a general method to write functional dilogarithm identities associated 
to the TBA formulation of the renormalization group evolution of integrable quantum 
models. It is worth to observe that such identities define quantities of topological na- 
ture, because they do not vary under arbitrary deformations of the vertex function 
j —I- Xj, provided that the initial ordering of the set {xj} is preserved. Our geomet- 
rical description of these identities in terms of triangulations (or quasi-triangulations) 
of three-dimensional manifolds provides us with the general solution of the recursion 

10 



relations associated to the Y systems and hence a simple explanation of their period- 
icity. Another consequence of this approach is to show that the y's can be expressed 
in terms of a unique periodic (or quasi-periodic) function x{6), related to Xj through 
Xj = x{9 + in-), where g depends on the model. In particular, it is easy to find the 
x{9) corresponding to the stationary (i.e. 9- independent) version of of the Y systems 
and their associated dilog identities. For instance one can verify, using eq.s ( ^.25 , 2^26| ) 



on iv+l 

that the stationary solution of the An system is given by x{9) = e ^+'^ and that of 
the D^ system is h{9) = 9. From a physical point of view, these functional identities 
define a quantity which is conserved during the renormalization group evolution of the 
system. In particular in the UV limit the system is described by a conformal field the- 
ory, then the r.h.s. of such identities must be expressible only in terms of conformal 
quantities such as the effective central charge and the conformal dimension A of the 
perturbing operators. For instance, in the purely elastic ADE related models the integer 



n of eq. (2.12) is given by the ratio -^- The RG flow of any quantum field theory can 



be viewed as a typical dissipative process, because the integration of degrees of freedom 
responsible of the trajectory between the UV and the IR fixed points yields an infor- 
mation loss of the system. Then it is amazing that there are conserved quantities along 
these trajectories. It would be very interesting, also in connection with the c theorem 
of A.B.Zamolodchikov, to find a field-theoretic description of such an invariant. 
Acknowledgements: RT thanks Michel Bauer for helpful discussions. 
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